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Background 


This is a report on the second six-month period of this NASA Grant. The 
purpose of the research is to study the integration of structural, dynamic, 
and aerodynamic considerations in the design-optimization process for helicop- 
ter rotor blades. This is to be done in three phases. The first phase, 
which represents the first year of effort and which is now completed, concen- 
trates on improving the structural optimization of our previous work to make 
it more realistic and to make it more compatible with the integration of 
dynamic and aerodynamic effects. A second goal of the first year's effort is 
to make preparations for years two and three through the development of a 
performance and loads program. 

Technical Personnel 


Four technical people are working on this project. Three are being paid 
through the grant, and one is being supported through institute funds. The 
principal investigator, David A. Peters, has contributed 0.9 man-months over 
this six-month period (15%). A Ph.D. student, Mnaouar Chouchane, has contrib- 
uted 2.0 man-months; and a post-doctoral fellow, Y.P. Cheng, has contributed 
0.48 man-months. In addition, an M.S. student, Mark Fulton, has worked on 
the project since September 15 while being supported by a Georgia Tech 
Research Assistantship. 

Scope of Effort 


The statement of work in our original proposal lists three tasks for the 
first-year effort (Phase I) with the third task being the largest. Task I 
is to bring on-line computer codes that could perform the finite-element 
frequency analysis of rotor blades. This has been done, as we summarized in 
our first status report (April 1987). Rather than bring in an existing code 
that might contain "extra baggage" not needed for our effort. Dr. Cheng wrote 
his own code specifically geared to our purposes. Table 1 summarizes the 
major features of this program. Although the program is based on the theory 
of our previous work, this new program is advanced in several ways, as 
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outlined in the Table. Of particular importance is the ability to model 
blade discontinuities and exact hinges and links, not previously available. 
This program is fully documented and operational. Thus, Task I is complete. 

The second task in the statement of work is to bring on-line an optimiz- 
ation code for the work. We have tried several (including OPT and CONMIN) 
and we have decided to use CONMIN. One great advantage of using this program 
is that it was also used in our previous work. Thus, we have a firm founda- 
tion of experience with it. Presently, we are using CONMIN routinely in all 
of our optimization studies. Thus, Task II is also finished. 

This brings us to the heart of our research, Task III, which involves 
the reproduction, refinement, and extension of our previous work. The 
original proposal lists four goals for Task III. The first is to add expli- 
cit volume constraints on the thicknesses and lumped masses used in the 
optimization. This was not explicitly carried out in earlier work. The 
second goal is to apply the specific aeroelastic constraint that the center 
of mass must be forward of the quarter chord in order to prevent flutter. A 
third goal is to include the bending-torsion coupling due to cg-ea offset 
within the blade cross-section. The fourth goal is to include some very 
simple stress constraints. In particular, the proposal calls for the axial- 
stress constraint and for a simple static bending constraint. (Vibratory 
stress constraints are reserved for the second and third years.) We have 
addressed all of these goals, have completed the first three, and are in the 
midst of completing the fourth. In addition, as a result of trying to meet 
the fourth goal, we have identified an additional area of research that must 
be completed before we can complete the fourth goal. The discussion below 
details the work done including this new task. 

Details of Added Optimization Task 

We begin with the new research task that is crucial to this work. This 
task is the use of optimization to match blade section geometry (thicknesses, 
material properties, etc.) to the structural beam properties of each cross- 
section (El's and GJ's). In our earlier work with NASA, the thickness of the 
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box beam at each section was chosen somewhat arbitrarily based on drawings of 
the blade. Then, any discrepancy between the beam properties of the box beam 
and those of the true blade were adjusted by the addition of lumped stiff- 
nesses and masses to make up the difference. The blade optimization was done 
on the box-beam properties, and the lumped "correction" properties were 
assumed to be unchangeable. An exception was that torsional stiffness for 
which the equivalent thickness of a second cell had to be determined by trial 
and error. This process consumed a great deal of man power, and we decided 
in our new research effort that it should be automated. Thus, Dr. Cheng has 
developed an optimization package that chooses cross-sectional properties in 
such a way as to match manufacturer data on blade beam properties. This not 
only allows beam properties to be converted into physical properties, but it 
also allows the possibility of optimization based on beam properties rather 
than on cross-sectional dimensions. In other words, one could optimize with 
El's and GJ's and then convert this optimum design to thicknesses; or one 
could optimize with geometrical properties as design variables. 

Table 2 outlines the present manner in which this process takes place. 
First, we try to minimize an objective function based on the differences 
between desired and actual beam stiffnesses. At each step, constraints are 
placed on stiffness wandering (x.); and these constraints are tightened as we 
approach the desired values. In step 2, we choose the shear modulus of the 
skin in order to match torsional stiffness. Finally, we adjust the lumped 
masses and the skin density in order to optimize a cost functional based on 
blade inertial properties. However, as we will see in results later in this 
report, this procedure has a weak point. In particular, although we can 
easily match stiffness properties, the method occasionally puts too much 
material into the blade skin and not enough into the box beam. Although this 
poses no particular problem for the frequency constraints, it does pose a 
problem when we add stress constraints. When the box-beam thickness is 
traded off for skin thickness, then the centrifugal stresses become too large 
in the box beam, since the skin does not carry this load. This has hampered 
the application of stress constraints because it makes the initial design 
unfeasible. 
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As a result, we are now experimenting with a revised procedure, outlined 
in Table 3. In this procedure, the skin properties are considered uniform 
along the blade. The thicknesses of the box beam are then used to match the 
four beam thickness properties. Because there are three thicknesses and four 
properties, we probably will need to include the trai ling-edge strip as an- 
other variable. Step 2 of Table 3 indicates that there may be a need to 
revise skin properties if the result is unacceptable, but this remains to be 
seen. Finally, we need to choose lumped masses to match inertia properties. 
In our new formulation, we have also recognized that rotary inertia has 
little impact on bending. Thus, we will only try to match the total polar 
moment of inertia rather than 1^ and I separately. Again, in the inertial 
portion we only have two free variables that must be used to match three 
properties. Therefore, we may need to add one additional parameter. 

Completed Optimization 


Although we have had some problems with the "property optimizer" men- 
tioned above, we have been able to complete most of the goals in Task III. 
This has been done on the Hughes blade since it is the most common type of 
articulated rotor. Table 4 shows the various stages used in this optimiza- 
tion process. In each stage, we include the volumetric constraint on cross- 
sectional properties (Fig. 2), we include the aeroelastic constraint on 
center of mass, and we include the torsional coupling due to cg-ae offset. 
The first page of this table illustrates that we begin with frequency-place- 
ment as the cost functional and bring the 3rd flap, 2nd inplane, and 4th flap 
into their desired windows. (Note that the 1st and 2nd flap and 1st inplane 
and torsion are already in their windows.) In the second phase, we bring the 
5th flap and 2nd torsion into their windows. With all frequencies now within 
the constraints, we have a feasible solution to begin the weight optimization 
(third page of Table 4). 

Table 5 shows the results of this three-phase optimization when the 
structural beam properties (El, GJ, etc.) are used as design variables. Note 
that the constraint on autorotational inertia prevents any great weight 
savings after the optimization. This is to be expected because the goal of 
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the optimization is to place frequencies. The minimum-weight requirement 
simply prevents the optimizer from doing this by adding unnecessary mass. 
Table 6 shows the same optimization when the physical properties (thickness- 
es, etc.) are used as design variables. One can see that there is little 
difference between the two methodologies except that the 6th flapping fre- 
quency is moved to a different window. The second method does require fewer 
iterations. Figure 2 shows the final physical properties from the two 
methods, and Fig. 3 shows the final structural properties. One can see that 
the optimization based on thicknesses gives more radial variation in proper- 
ties. We believe that this is due to the problems mentioned earlier with 
respect to property selection. One can also see the large variation in 
thicknesses around the r=50 station. This is also due to this selection 
problem. It is the very small thicknesses that give difficulties when we 
attempt to apply stress constraints. We believe that these problems will be 
eliminated with our new strategy, and we expect to be able to add the stress 
constraints shortly. 

It should also be noted here that the inclusion of bending-torsion 
coupling caused other problems that we had not anticipated. In particular, 
because of frequency cross-overs, it was sometimes hard to track a particular 
mode shape and frequency to see which frequency was really in which window. 
This is one reason that we separated the frequency placement into three sub- 
problems (rather than only two as in our previous work). This helped us to 
track modes more efficiently. 

Aerodynamic Loads and Performance 

In our first status report, we described the development of a perform- 
ance and loads program that would include stall and that could be used in the 
second and third years of our study. We have continued work on this program 
with Mr. Chouchane taking the lead. In our efforts to push the program to 
higher advance ratios (and higher loads), we discovered that we could not 
trim the helicopter beyond an advance ratio of 0.30 to 0.35. In order to 
find the reason for this limitation, we re-examined the optimization that 
previously produced gains and time constants for the trim algorithm. Based 


5 


on this, we found that the previous solution, although adequate at low 
speeds, did not give optimum performance when high speeds and/or stall were 
present. Therefore, we have been attempting a new optimization. Figures 
4a-4c show time histories of rotor controls at an advance ratio of 0.25. 
Five different combinations of gains (K0 and Kl) and time constants (TO and 
Tl) were used. The objective is to reach a steady-state as soon as possible 
and with as little steady oscillation as possible (<+0.3°). Based on this, 
the "+" curve seemed the best of the five. Figures 5a-5c then provide the 
effect of changing the gain Kl alone. Based on this, Kl=.04 seemed best. 
Consequently, we have gone ahead with these values. However, another stu- 
dent, Mark Fulton, has been given the task of doing a formal optimization to 
find the absolute best gains for this system. He is supported by Institute 
funds. 

Once these sub-optimum gains were determined, we were able to locate the 
source of our trim problems. The rotor coupling matrix in our trim algorithm 
was simply not accurate enough at high advance ratio. Thus, we replaced it 
with a more accurate set of formulas based on the work of Peters and Ormis- 
ton. Figures 6a-6c show the results at advance ratios from 0.25 to 0.39. 
The new coupling matrix has pushed results up to y=0.375. Beyond that, even 
these advanced formulas are inadequate. Figures 7a-7b show blade response 
for these same cases. They show that a periodic solution with no 1/rev has 
indeed been obtained up to y=0.375 but not beyond. Therefore, we are develop- 
ing an adaptive control that will optimize itself during the trimming process 
to match the particular couplings of any flight condition. This should allow 
us to develop the loads data we need in future phases. 

Future Work 


Our primary objective now is to complete the the new cross-section 
optimizer so that we can apply the stress constraints. As stated in our 
proposal, these initial stress constraints will be based on static axial 
loads and on the static lift necessary to lift the aircraft (flap bending 
only). Our second objective is the completion of Phase II which includes the 
addition of performance constraints within the optimization process. In this 
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case, the weight objective function will be replaced by a combination of 
weight and performance which will result in the best payload capability for a 
given power. Although we will begin with hover, we cannot stop there, 
because most helicopters are not optimized for hover performance alone. 
Therefore, we will need to include forward-flight performance. Our trim and 
stall methodologies, developed in Phase I will be useful in this context. 

Finally in Phase III, we add vibratory loads and vibratory stress con- 
straints. Here, again, we plan to relate everything to a performance index. 
In particular, higher stresses imply stronger ( i . e . , heavier) blades; and 
higher vibrations also imply vibration suppression devices which are a weight 
penalty, as well. As for the source of these vibratory loads, we will begin 
with a simple load spectrum. Towards the end, however, we will apply the 
loads program on which we are now working. 
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Table 1: Some Features of the Finite Element Program of Rotor Blade. 
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Gradient information is calculated analytically. 

. Subspace iteration method is used to increase the efficiency of the iteration process. 



Table 2: Procedures Now Used to Determine Section Properties from Given Structural Data. 
* Step 1: Design variables: t. s p a v or piskin thickness) 
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( i-Xg)(e 0 )< e <(l+t 7 )(e 0 ) 

Where wj is weighting scalar, if is tolerant value, ( )o is given data, ( )f represents quantity i 
flapping direction, ( ) c represents quantity in chordwise direction, and e is elastic offset. 



Table 3: Proposed Procedures to Determine Section Properties from Given Structural Data. 
* Step 1: Design variables: t, s v s 2 ( Given skin thickness & G(skin) ) 
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( 1-X 4 ) ( plfo X plf <( 1 +t 3 ) ( plfo ) 
( 1-Xg ) ( pic o)< pic < ( 1 + tg ) ( pl c 0 ) 
d-XgHeoX e <(I+T 7 )(e 0 ) 



Table 4: Procedures of Optimizing the Hughes Articulated Rotor Blade. 
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(blade natural frequency) / (rotor rotational frequency). 



Table 4: Procedures of Optimizing the Hughes Articulated Rotor Blade. (Continue) 

BaSed 0n: Structural Properties Physical Properties 
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Objective ( [ p ( flapping-5th) - 11.5 ] 2 

+ [ p( torsion-2nd ) - 12.5 ] 2 ) 



Table 4: Procedures of Optimizing the Hughes Articulated Rotor Blade. (Continue) 
Based on Structural Properties Physical 
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constraint ( elastic offset > -0.54 ). 



Table 5: Optimization Results of Hughes Blade ( Structural Properties as Design Variables). 
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Iterations(CONMIN) 
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Iterations(CONMIN) 
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Figure 1: Typical Blade Section 
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Material properties: Structural data: 

Density (box beam, lumped weights, & honeycomb filler) EA, Elf, EIc, GJ 

E & G (box beam & skin) D If. D L e 








Figure 2: Structural Data Comparison (Continue) 
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Figure 3: Section Thickness Compari 
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Figure 3: Section Thickness Comparison (Continue). 
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figure 4a: AUTO-PILOT CONTROL MU = .25 
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figure 4b: AUTO— PILOT CONTROL MU = .25 
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FIGURE 4c: AUTO — PILOT 
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figure 5a: AUTO-PILOT CONTROL 

MU = .25 T0=3. 141 6 T1=9.425 
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figure 5c: AUTO-PILOT CONTROL 

MU = .25 T0=3. 1416 T1 =9.425 



40.00 60.00 80.00 100.00 120.00 140.00 

AZIMTH ANGLE 









figure 6b: AUTO— PILOT CONTROL 
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figure 6c: AUTO-PILOT CONTROL 
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FIGURE 7b: FLAPPING RESPONSE 
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